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A study of the bound states for square potential wells with position-dependent mass 
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A square potential well with position-dependent mass is studied for bound states. Applying 
appropriate matching conditions, a transcendental equation is derived for the energy eigenvalues. 
Numerical results are presented graphically and the variation of the energy of the bound states are 
calculated as a function of the well-width and mass. 
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I. INTRODUCTION 

The concept of position-dependent mass comes from 
effective-mass approximation of many-body problem in 
condensed matter physics [H 0, 0, H IE lEM 0- In re- 
cent times a good number of articles have been published 
[ISIIllIlHQEEIIlElliin this field. The 
Schrodinger equation with position-dependent mass has 
been studied in the contexts of supersymmetry, shape- 
invariance, Lie algebra, point-canonical transformation, 
etc. It is well known that the kinetic energy operator in 
this case belongs to the two-parameter family yj 

T(x) = — (m a pm^pm y + m 7 pmP pm a ), (Id) 

where m = m(x) and p — —ihd/dx, with the constraint 



a + /3 + 7= -1. 



(1.2) 



However, the correct values of the parameters a, 8, 7 for 

a specific model is a long-standing debate 

For example, in the case of a potential step and barrier 

, and in some one-dimensional potential wells [j| with 

varying mass, the following kinetic energy operator was 

chosen 



T ( x ) = \\ P — P 
2 \ m 



(1.3) 



which corresponds toa = 7 = 0, j3 = — 1. Then, it was 
shown that the reflection and transmission coefficients 
Q , as well as the discrete spectrum , are different com- 
pared to the constant-mass problem. 

In this article we are going to study the following phys- 
ical problem: the potential energy has the form of a well, 



both in symmetric as well as asymmetric form, and the 
kinetic energy is given by Hl.ljl with a mass function m(x) 
which has different constant values inside and outside the 
well. This problem has also interest from the point of 
view of applications in carbon nanotubes and quantum 
dots, as can be seen in [2(],E1- Our purpose is first to 
find an ordering of the kinetic energy term appropriate to 
this problem, and then, to study in full detail the bound 
states of this model and to compare the results with the 
conventional constant-mass case. 

Several problems which resemble our model have been 
considered previously in the literature. For example, the 
scatterin g st ates of a potential step or barrier were stud- 
ied in [sl \22\ . Other special square well potentials were 
analysed in [23], and the bound states for some finite 
and infinite wells were studied using other kinetic en- 
ergy operator and matching conditions in Q. However, 
our approach is different from the very beginning be- 
cause we propose an ordering for kinetic term based on 
a specific argument for the matching conditions. This 
ordering is also used in |24| for studying connection rules 
for effective-mass wave functions across an abrupt het- 
erojunction. 

The structure of this article is as follows. Sec. [H] is 
devoted to fix the ordering in the kinetic term. In Sec. lIIII 
a transcendental equation determining the energy values 
for the bound states in the case of a potential energy well 
is derived. The numerical results are shown graphically 
and discussed in Sec. II VI Finally, Section Ivl contains the 
conclusions of our work. 



II. MATCHING CONDITIONS AT THE 
DISCONTINUITIES OF THE MASS AND 
POTENTIAL FUNCTION 
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The Hamiltonian operator for the position-dependent 
mass problem is given by 

H = T{x) + V{x) (2.1) 

where T(x) is the operator defined in Ijl.lfl and V(x) 
denotes the potential term. The one-dimensional time- 
independent Schrodinger equation for the stationary 



2 



states is 



and Eq. I|2.6[) takes the following simple form 



Hil){x) = Eil){x). 



(2.2) 



Substituting and (|2.1|) into (|2.2|) and taking units 

such that ti 2 — 2, we get the Schrodinger equation for 
the generalized kinetic energy operator 



d 2 ip to' dip 
dx 2 m dx 



1 / m to' , 
- [v v — +m(E-V) 

2 \ to to/ 



where 



r] = a(-f + 2) — 7(a + 2), ^ = a + 7, 



= 

(2.3) 
(2.4) 



and m' = dm/dx. If in Eq. (|2.3(l we make the following 
transformation 



ip(x) = m{x) 1 ^ 2 <p{x) 



(2.5) 



we can eliminate the first order derivative of ip with re- 
spect to x, to arrive at 



d 2 ^ 



+ 



2to 



Tj \ TTt 



1 2 



+ TO (£ - V) 



0=0. 

Our aim is to generalize the usual matching condi- 
tions for the wave function ip(x) suitable for this position- 
dependent mass problem. We will focus on the case of 
mass discontinuities. 



A. Matching conditions 

Now, let us assume that the mass m(x) has a finite 
discontinuity at x = a of the form 

m{x) ~ mi(x)Q(~x + a) + m2(x)<d(x — a) (2.7) 

where mi(x) and to-2(x) are smooth functions and the 
unit step function <d(x) is defined as 



0(3!) 



1, x > 
0, x < 0. 



(2.8) 



If we use <|2.7[) in Eq. I|2.6|) , we see that it leads to strong 
discontinuities at x = a as well as some terms that re- 
quire a careful interpretation. For instance, we should 
take into account that 0' (x) = d<d(x)/dx = S(x), where 
S(x) denotes the Dirac delta distribution. In order to 
eliminate these problems in Eq. I|2.6I) . we can choose 



0. 



0, 



(2.9) 



where 77 and v were given in (|2.4|l . From these conditions, 
we obtain the values of the parameters a — 7 = — 5 and 
(3 = 0. With such values, the kinetic energy operator 
(11.11) becomes 



T 1 ( 1 2 

2 U/m 



1 



(2.10) 



-^ + m {V~E)cp 



(2.11) 



Now, it is easy to get the matching conditions for ip at 
x = a. To do this, first we integrate Eq. I|2.11|l around 
the discontinuity point x — a 



<j>'(a + h)-<j)'(a-h) 



a+h 



m(x) (V(x) — E)(f>(x) dx. 



-h 



(2.12) 

In the interval (a — h, a+h), the functions m(x) and V(x) 
have finite discontinuities at x = a and 4>(x) is bounded. 
Therefore, when h goes to zero, the integral at the r.h.s. 
of 12.12fl tends to zero. This means that <fi'(x) (and also 
4>{x)) will be continuous at x = a. In conclusion, using 
the transformation (|2.5(l we have arrived at the following 
matching conditions: the wave function ip{x) is such that 
the two functions ip(x) / •*/ m(x) and (ip(x) / -\J m(x))' are 
continuous at x = a. Mathematically it may be expressed 



as 



ip(x) 



y/m(x) 



ip(x) 



x—a—0 



4>{x) 
V m ( x ) 



\/m(x) 



tp(x) 



x—a-\-0 



ijj(x) \ 



=a+Q 



y/m(xj 



\Jm{x) t 

(2.13) 

where a is an interior point in the domain of the problem. 
It is straightforward to check that the above conditions 
are consistent with the definition of a time-independent 
inner product in the space of square integrable eigenfunc- 
tions (ip(x),<f>(x)) — J ip(x)*tfi(x) dx, as well as with the 
conservation of the current density 



j(x) = -i ip(x) 



1 dip(x) dip(xy 



m{x) dx dx m(x) 

(2.14) 

It is clear that these matching conditions are consistent 
with the Hermitian character of the Hamiltonian 12.1(1 
with (|1.1J) . We should stress that the above mentioned 
boundary conditions indeed define a self-adjoint problem 
(this issue will be addressed in detail elsewhere j25|). 



III. A SQUARE POTENTIAL WELL AND A 
STEP MASS 

A. Asymmetric well 

Let us consider an asymmetric well of the form 



Vi, x < —a 
n.r) = { 0, \x\ < a 
V2, a < x 



(3.1) 



3 



with the position-dependent mass 

mi, |x| > a 
m 2 , \x\ < a 



m(x) 



(3.2) 



where mi, m%, V\, and V 2 are constants such that V% > 
Vi, mi,TO2 > 0, and mi ^ m 2 . Next, we will study 
the bound states (0 < E < Vi < V 2 ) for this problem 
using to the matching condition obtained in the previous 
section. 

The Schrodinger equation (|2.3|l for the wavefunction 
ip(x) has the following form in each region 

d 2 ip 



— |- - k\i> = 0, fci = y/rm^-E), x < -a, (3.3) 



d 2 ip 
dx 2 



k 2 ip = 0, k 2 = \Jm 2 E, 



\x\ < a, (3.4) 



-^-^ = 0, k 3 = y/rrnty-E), x > a. (3.5) 
The physical solutions of these equations take the form 



Ae klX , 



x < —a 



ij){x) = { Csm(k 2 x + 8), \x\ < 



(3.6) 



Be'" 3 



k 3 :r 



x > a 



where A, B, C, and 9 are constants to be determined by 
the boundary and matching conditions. Using the con- 
tinuity conditions (|2.13|) for this solution at x = —a, we 
obtain the following two equations 



A 



m 2 
mi 



1/2 



- ki a 



= Csin(-k 2 a + 9) (3.7) 



/ \ 1/2 

A \nVj kie ~ kia = Ck 2 cos(-k 2 a + 9) (3.8) 

and at x = a, we get 
/ \ 1/2 

B \m 1 ) C ~ k3a = Csin ( fc 2a + 60 (3.9) 
/ \ V 2 

B[—\ k 3 e~ k3a = Ck 2 cos(k 2 a + 8). (3.10) 
\mij 

From these four equations, we find 

ki — k 2 cot(— k 2 a + 9), fc 3 = — k 2 cot(k 2 a + 9) (3.11) 
or 

k 2 



s\n(—k 2 a + 9) 
sin(fc 2 a + 9) 



miVi-k 2 {^-l) 



HV 2 -k 2 (^-l) 



(3.12) 
(3.13) 



Now, if we eliminate 9 from l|3.12|l and (|3.13(l , we get the 

transcendental equation 

k 2 



2k 2 a = nir — sin 



miV 2 -k 2 (^-l) 



sm 



sJmiVi-k 2 {^-l) 



, (3.14) 



where inside the well, the momentum k 2 of the bound 
states is obtained from the values n ~ 1,2,3,..., and 
the range of the inverse sine is taken between and ^ . 
Since E = (k 2 ) 2 /m 2 , the roots of this equation also give 
the energy values of the bound states. 

As the maximum value that k 2 could take is \J m 2 V\ , 
from Eq. (|3.14|l . we get an inequality for the number of 
bound states 



2a \j rn 2 Vi > — — ^ tt — sir 



\Jm 2 Vi 



\/m 2 Vi + miAV ' 

(3.15) 

with AV = V 2 - V x > 0. This means that the total 
number of bound states TV will be the highest n satisfying 
this inequality. When we consider (|3.15|) as an equation 



for TV = 1,2,.. 



2a^/m 2 Vi 



N 



V m 2 Vi 



\Jm 2 Vi + miAV 

(3.16) 

we say that the parameters of the well are "critical" , in 
the sense that by modifying slightly their values we will 
have one bound state less or one bound state more. 

From (|3.16|) . for N = 1 we see that the first bound 
state of the asymmetric square well will appear when the 
following condition is satisfied: 



sin 



\Jm 2 V~i 
\Jm 2 Vi + miAV 



= J - 2ayjm 2 Vi . (3.17) 



Therefore, we can say that for V\ ^ V 2 , there are al- 
ways some values of the parameters of the well that do 
not allow for bound states, as it is also the case for the 
conventional constant-mass problem. 

The influence of mi (the mass outside the well) on the 
number of bound states is not so important. Let us call 
/(mi) the function on the l.h.s. of Ij3.17|l . then 



lim /(mi) 

mi— >0 



lim /(mi) = - 

mi — *oo Z 



(3.18) 



This means that as we increase mi the spectrum of bound 
states either will remain the same or will have one value 
less. This is in sharp contradistinction to the influence of 
m 2 , that can change any number of levels in the discrete 
spectrum. 

The conventional constant-mass case is obtained by 
taking mi = m 2 — m in Eq. lj3~Hjl to get the well 
known formula 26] 



2k 2 a 



nn — sm 



\JmV 2 



\JmVi 



(3.19) 



where n = 1,2,3,... The inequality for the number of 
bound states is obtained now from l|3.15|l if we put mi = 

m 2 = m: 



2a\JmVi > 



n I ;T — sm 

2 



— -. (3.20) 
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However, we must stress that in the position-dependent 
mass case the critical values given by formula 13.15( 1 for 
the number of bound states depend on both mi and 7712, 
in particular, as m 2 — * 0, we see from (|3.17() that no 
bound states will remain in the well. 

B. Symmetric well 

If Vi = V 2 = V, from (|3.11l) . we get the energy equa- 
tions corresponding to even and odd eigenfunctions by 
replacing 9 = tt/2 and 9 — 0, respectively 

— !- = tanfc2a, 7— = — tan k 2 a. (3-21) 

k 2 ki 

In this case, we can also write for both even and odd so- 
lutions the transcendental Eq. 13. 141 1. using k 2 = ^/rn^E, 
in the following way 

9fc 2 » = mr- 9 sin' 1 — k2 , (3.22) 

^kj + mi (V-E)' 

where n = 1,2,3,... Here, the argument of the inverse 
sine is always less than 1 for V — E > 0, and therefore 
the formula is well defined. If we make E = V in (|3.22(l 
the number of bound states N is given by the highest n 
satisfying the inequality 

2a^m 2 V > (n— 1)tt. (3.23) 

This means that there are always bound states for the 
symmetric well. Even more, once the depth V and the 
width 2a of the well are fixed, the number of bound states 
does not depend on mi , it depends only on the mass in- 
side of the well (7712). The relation (|3.23(l for the number 
of bound states coincides with that obtained also by Plas- 
tino et al using different matching conditions. However, 
the energy equations (|3.21() - (|3.22l) are different from 0- 
In fact, when the inside mass is bigger (smaller) than the 
outside mass, then our energy spectrum is lower (upper) 
than the spectrum given in |9(. 

The critical values of the well (which determine when 
we have one bound state more) are obtained from 

2ay/m 2 V = (N-l)n. (3.24) 

Thus the critical values of a are linear in the number 
TV — 1 of bound states, while the critical values m 2 grow 
as the square of (TV — 1). In conclusion, the important 
parameter here is the mass m 2 inside the well, while the 
mass outside (mi) plays a secondary role. 

The number of bound states for the conventional 
constant-mass symmetric well are given also by (I3.23f) 
replacing m 2 — > m. However, the energy of these bound 
states is defined through 1(3.22(1 by choosing mi = m 2 — 
m 

k 2 a = n sin -1 — 2 (3.25) 

2 y/mV 

Hence the values of the energies will be different from the 
position-dependent mass. 



IV. NUMERICAL RESULTS AND DISCUSSION 

In this section we will discuss some numerical results 
and graphics obtained from the formulas of the previous 
sections, paying attention to the consequences of mass 
variation. 

Consider the situation of a particle of mass mi in a 
square well potential, where we have altered the condi- 
tions inside the well to produce a different effective mass 
m 2 . We want to study the effects that this change of 
mass will produce on the bound states. 

A. Features of the asymmetric well 

The main characteristics of the asymmetric well may 
be summarized as follows: 

1. If Vi < V 2 , it can be seen from formula (|3.15|) that 
the number of bound states depends linearly on the 
width 2a of the asymmetric well as shown in Fig.^ 




1 2 3 4 5 6 



a 

FIG. 1: Plot of the energy values and the number of the bound 
states depending on the width of the well 2a when the other 
parameters are fixed as: Vi = 1, V2 = 2, mi = 1, U12 = 2 
(solid lines), ni2 = 1 (dotted lines), and mi = 0.5 (dashed 
lines). 

where we have plotted three situations: 

(a) the usual constant-mass case, mi = m 2 = 1, 

(b) mi = 1, m 2 = 2, that is the case where the 
mass inside is bigger than outside the well, 

(c) mi = 1, m 2 = 0.5, that corresponds to a 
lighter mass inside the well. 

Thus, the effect of m 2 on the bound states is quite 
strong for any value of a, as can be seen in Fig. ^ 
When m 2 is decreasing, the number of bound states 
also decrease, while each energy level is higher. In 
Table it is shown the effect of these three values 
of m 2 on the first critical values of the width 2a. 

2. In Fig. |21 we plotted the bound state energies as we 
change only the inside mass m 2 for mi = 1 and 
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TABLE I: The critical values of a obtained from Eq. 13.161 . 
determine the width where a new bound state appears in the 
well. 



mi 


m2 


a,™ 


a< 2 > 


a< 3 ) 


a< 4 > 


a( 5 > 


a< 6 > 


1 


2 


0.2176 


1.3283 


2.4390 


3.5498 


4.6605 


5.7712 


1 


1 


0.3927 


1.9635 


3.5343 


5.1051 


6.6759 


8.2469 


1 


0.5 


0.6755 


2.8970 


5.1184 


7.3398 


9.5613 


11.7827 




10 20 30 40 50 60 

ni2 



FIG. 2: Plot of the energy values and the number of the bound 
states depending on the mass m,2 when the other parameters 
are fixed as: a = 1, Vi = 1, Vi = 2, mi = 1 (solid lines), and 
mi = 10 (dashed lines). 

mi = 10 (leaving the other parameter fixed). We 
see that the number of bound states is almost pro- 
portional to the square root of 7772 due to l|3.15[l . 
while the energy value of each bound state de- 
creases with 777-2 ■ The influence of V\ is similar as 

that Of 777,2- 

3. In the same figure, we can see that the influence of 
777 1 (the mass outside the well) on the number of 
bound states is not so important. Table ILTI shows 
the influence of mi on the first critical energy values 
of 7772- The parameter V% has similar qualitative 
effects on the energy values and the number of the 
bound states as the parameter m\. 



B. Features of the symmetric well 

In the symmetric well the importance of the mass 777,2 
inside the well is even stronger than in the asymmetric 
well. The main characteristics of the symmetric well arc 
the following: 

1. The number of bound states depends linearly on a, 
the square root of the 777,2 and V, as in the asym- 
metric well, but it is independent of mi. 

2. The dependence on 7772 is shown in Fig. |3 a detail 
is given in Fig. As a consequence, the critical 



TABLE II: The critical values of m,2 obtained from Eq. (13.161 , 
determine the mass where a new bound state appears in the 
well. 



mi 




(2) 

m 2 


(3) 


(4) 


(5) 


(6) 


1 


0.2899 


3.3189 


10.8186 


23.1821 


40.4642 


62.6758 


10 


0.4618 


4.2715 


12.3087 


24.9461 


42.3710 


64.6628 




5 10 15 20 25 30 35 

/77 2 



FIG. 3: Plot of the energy values and the number of the bound 
states depending on the mass wi2 when the other parameters 
are fixed as: a = 1, Vx = V2 = 2, mi = 1 (solid lines), 
mi = 10 (dashed lines), and mi = rri2 = m (dotted lines). 
The critical values for the three cases are the same: m^ — 0| 
m 2 2) = 1.23, m 2 3) = 4.93, m 2 4) = 11.10, m 2 5> = 19.74, ■m { 2 ) = 
30.84, obtained from Eq. fTZR . 



values of 7712 are the same, independently of the 
values of m\. 

3. The energy values given by (|3.22[) are only slightly 
affected by the variation of mi, the mass outside 
the well. This is shown in Fig. [3| for the values 
mi = 1, m\ = 10 and toi = 7772 = m. A detail of 
this plot is shown in Fig.0J where the energy levels 
of the case mi = 1 coincide with the conventional 
constant-mass case (mi = 7772 = m) at the point 
7772 = 1- Clearly, the energy levels of the case mi = 
10 will also coincide with those for constant-mass 
case mi = 7772 = m at the point 7772 = 10 (see 

Fig. El). 



V. CONCLUSIONS 

In this article we have studied a potential well both 
in symmetric and asymmetric form, with different con- 
stant mass inside and outside the well. We have proposed 
specific values for the ordering parameters in the kinetic 
term based on the simple argument that strong singular- 
ities in the effective-mass Schrodinger equation should be 
avoided. However, this choice is not unique: for instance, 
the symmetric well has been considered in |9| with a ki- 
netic term proposed in Q. Similar to the conventional 
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FIG. 4: Plot of the energy values and the number of the bound 
states depending on the mass mi when the other parameters 
are fixed as: a = 1, V\ = V% = 2, mi = 1 (solid lines), 
mi = 10 (dashed lines), and mi = m-2 = m (dotted lines). 

constant-mass case we have obtained a transcendental 
equation for bound state energies. We have shown our 
numerical results for different values of the parameters. 
Many interesting differences from conventional constant- 
mass situation have been pointed out. In particular, the 
bound states are controlled mainly by the mass inside 
the well, while the mass outside, in general acts simply 
as a tuning of the specific values of the energies in the 
spectrum. The experiments have also shown that the en- 



ergy level spacing is very sensitive to the effective mass 
inside the well, while the effects of the other parameters 
are not so important j2£j • Finally, we have remarked the 
influence of the matching conditions, related to the ki- 
netic term, on the discrete spectrum by comparing our 
results with those of |9j ■ 

We have also determined some critical values of the 
well, that is the values of the potential and mass func- 
tions giving rise to the new bound states in the well. 
These critical values are important in the study of the 
scattering states, because they fix the conditions where 
the transmission coefficients reach the maximum values. 
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